arXivil 504.02368vl [quant-ph] 9 Apr 2015 


Optical hyperpolarization of 13 C nuclear spins in nanodiamond ensembles 


Q. Chen 1 ’ 2 , I. Schwarz 1 ’ 2 , F. Jelezko 2 ’ 3 , A. Retzker 4 and M.B. Plenio 1 ’ 2 
1 Institut fur Theoretische Physik, Albert-Einstein-Allee 11, Universitat Ulm, 89069 Ulm, Germany 
2 IQST, Albert-Einstein-Allee 11, Universitat Ulm, 89069 Ulm, Germany 
3 Institut fur Quantenoptik, Universitat Ulm, 89073 Ulm, Germany 
4 Racah Institute of Physics, The Hebrew University of Jerusalem, Jerusalem, 91904, Israel 

Dynamical nuclear polarization holds the key for orders of magnitude enhancements of nuclear 
magnetic resonance signals which, in turn, would enable a wide range of novel applications in 
biomedical sciences. However, current implementations of DNP require cryogenic temperatures 
and long times for achieving high polarization. Here we propose and analyse in detail protocols 
that can achieve rapid hyperpolarization of 13 C nuclear spins in randomly oriented ensembles of 
nanodiamonds at room temperature. Our protocols exploit a combination of optical polarization 
of electron spins in nitrogen-vacancy centers and the transfer of this polarization to 13 C nuclei 
by means of microwave control to overcome the severe challenges that are posed by the random 
orientation of the nanodiamonds and their nitrogen-vacancy centers. Specifically, these random 
orientations result in exceedingly large energy variations of the electron spin levels that render 
the polarization and coherent control of the nitrogen-vacancy center electron spins as well as the 
control of their coherent interaction with the surrounding 13 C nuclear spins highly inefficient. We 
address these challenges by a combination of an off-resonant microwave double resonance scheme in 
conjunction with a realisation of the integrated solid effect which, together with adiabatic rotations of 
external magnetic fields or rotations of nanodiamonds, leads to a protocol that achieves high levels of 
hyperpolarization of the entire nuclear-spin bath in a randomly oriented ensemble of nanodiamonds 
even at room temperature. This hyperpolarization together with the long nuclear spin polarization 
lifetimes in nanodiamonds and the relatively high density of 13 C nuclei has the potential to result in 
a major signal enhancement in 13 C nuclear magnetic resonance imaging and suggests functionalized 
and hyperpolarized nanodiamonds as a unique probe for molecular imaging both in vitro and in 
vivo. 


I. INTRODUCTION 

Nuclear magnetic resonance (NMR) GUI and mag¬ 
netic resonance imaging (MRI) Q have evolved to be 
powerful techniques to extract molecular-level informa¬ 
tion in a wide variety of physical, chemical and biological 
applications. Recently, 13 C based MRI has emerged as a 
new platform enabling essentially background-free imag¬ 
ing of non-proton nuclei with the additional possibility 
for the local and permanent destruction of the signal by 
means of radio frequency (RF) pulses [4]. However, the 
low signal-to-noise ratio for 13 C based MRI is not suffi¬ 
cient for most clinical and research applications due to 
the combination of a relatively low gyromagnetic ratio 
of 13 C and of the low natural abundance of this nucleus. 
Several strategies have been proposed to enhance the sen¬ 
sitivity of 13 C based MRI. Particularly promising in this 
context is the hyperpolarization of the 13 C nuclei, i.e. 
the generating of a large, non-thermal 13 C nuclear spin- 
polarization. Indeed, hyperpolarized 13 C based MRI has 
shown exciting potential for in vivo applications, espe¬ 
cially for novel metabolic imaging GHz|. One of the most 
powerful methods for the generation of hyperpolarization 
is dynamic nuclear polarization (DNP) (8l4l2|. in which 
a large polarization of electron spins is transferred to nu¬ 
clear spins, resulting in an enhancement of the MRI sig¬ 
nal by several orders of magnitude. 

Over the last decade, several breakthroughs have oc¬ 
curred in the manufacturing, surface treatment and use 
of nanoparticles for biomedical applications [Tj|. Bio¬ 


compatible nanoparticles present an attractive platform 
for hyperpolarization, as they can be functionalized for 
molecular specificity, ca exhibit long nuclear spin relax¬ 
ation times, and contain numerous nuclear spins of a sin¬ 
gle species (e.g. carbon). Long-lived hyperpolarization 
has been demonstrated in silica nanopartices 0 and 
nanodiamonds HE is], using standard DNP protocols, 
i.e. transferring the high thermal electron spin polariza¬ 
tion at cryogenic temperatures in a strong external mag¬ 
netic field to nuclei. Nanodiamonds in particular offer 
exciting possibilities as novel hyperpolarized probes. In 
addition to excellent biocompatibility 0 , the possibility 
for surface functionalisation and nuclear spin relaxation 
time of several minutes [l8| , nanodiamonds contain crys¬ 
tal point-defects with unique optical and magnetic prop¬ 
erties. Amongst these the negatively charged nitrogen- 
vacancy (NV) center stands out as its electron spin can be 
polarized within microseconds by optical pumping while 
exhibiting a relaxation time in the millisecond range even 
at room temperature [19]. This establishes the NV cen¬ 
ter electron spin as promising candidate for achieving 
nanodiamond hyperpolarization at ambient conditions - 
room-temperature optical DNP. 

Previous work has demonstrated that in bulk diamond 
electron spin polarization of the NV center can be gen¬ 
erated and subsequently be transferred via hyperfine in¬ 
teractions to nearby nuclei. This polarization can then 
be detected indirectly via the NV center [H, l20l427| or 
directly via an NMR scanner [HI- To this end, in bulk 
diamonds a high magnetic field is aligned with the quan- 
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tisation axis of the NV center to enable optical pump¬ 
ing to a specific electron spin state. Controlled inter¬ 
action of the electron spin with specific nearby nuclear 
spins species can then be achieved for example by the ap¬ 
plication of continuous wave microwave (MW) radiation 
applied to the NV center electron spin [24|. When the 
Rabi frequency of the microwave driving field matches 
the Larmor frequency of a specific nuclear spin species 
(achieving a Hartmann-Hahn (H-H) resonance [29]), flip- 
flops can occur between the microwave dressed states of 
the electron spin of the NV center and the surrounding 
13 C nuclear spins j30l - [32j as was recently demonstrated 
experimentally [24]. 

However, medical MRI applications require large en¬ 
sembles of nanodiamonds. In such ensembles (powder or 
solution) the angle between the external magnetic field 
and the natural orientations of NV centers and there¬ 
fore their quantisation axis is randomly distributed. This 
brings about several challenges for achieving high levels 
of hyperpolarization that need to be addressed by means 
of carefully designed protocols. First, the energy levels 
of the NV spin are then distributed over a large energy 
range of possible values so that the applied microwave 
field will be resonantly coupled to an exceedingly small 
fraction of the NV spins. Secondly, in an external mag¬ 
netic field, additional limitations concern the optical po¬ 
larization of the NV spins. In particular, the NV spins 
will be initialized to different states relative to the labo¬ 
ratory frame depending on their angle with the external 
magnetic field, thus resulting in a small net polarization. 

In this work we will show to address these challenges 
to achieve a radical enhancement of the hyperpolariza¬ 
tion in a nanodiamond ensemble. To this end we use an 
off-resonant driving and the integrated solid effect (ISE) 
for the NV center spin, resulting in a robust spin po¬ 
larization transfer between the NV and nuclear spins. A 
very large fraction of the NV electron spins can be coher¬ 
ently coupled to neighboring nuclear spins. Additionally, 
adiabatic rotation of the magnetic fields, or the Brown¬ 
ian rotation of the nanodiamonds themselves (in powder 
or solution respectively), is proposed for extending the 
polarization scheme to almost the entire nanodiamond 
ensemble. We also discuss the effect of spin diffusion in¬ 
duced by the nuclear dipole-diploe interactions which will 
further support the polarization of large volumes. 

The paper is organized as follows. In Sec. [IT] we present 
the coherent coupling and polarization of a single nuclear 
spin near a randomly oriented NV spin. This serves the 
introduction of the principal challenges in nanodiamond 
ensemble polarization, the discussion of the robust initial¬ 
ization and construction of suitable dressed states for the 
NV center electron spin and the realisation of near reso¬ 
nant coupling with the nuclear spin by using off-resonant 
driving and the integrated solid effect. In Sec. Hill we 
introduce mechanisms to extend our polarization proto¬ 
cols to the entire nanodiamond ensemble, both for nan¬ 
odiamond powder and for nanodiamonds in a solution. 
In Sec. EO we consider our protocol for multiple 13 C 


nuclear spins and include the effect of dipolar coupling 
among the nuclear spins and in particular the benefits of 
nuclear spin diffusion. In Sec. El we demonstrate the ap¬ 
plicability of our polarization schemes for another very 
broad range of NV center orientations. In Sec. [VI] we 
take into account the effect of depolarization processes 
on the efficiency of our polarization protocol. Finally, 
the discussion and conclusion parts are given in Sec. IVIII 
and Sec. EUD respectively. 


A Z 



FIG. 1. Schematic of dressed-state resonant coupling of 
NV spins and nuclear spins in a nanodiamond ensemble and 
energy-level diagrams, (a) The random orientations of the NV 
spins are uniformly distributed over the unit sphere. Small 
yellow circles represent examples with orientations indicated 
by red arrows which represent the unit vector pointing from 
the nitrogen to the vacancy of an NV center (see the example 
of such an NV center in the dashed circle, the axis pointing 
along ‘N-V’ forms the natural quantization axis), (b) Ground 
electronic spin states of an NV center in a strong external 
magnetic field. Microwave driving fields are applied to off- 
resonantly achieving an effective resonant double quantum 
transition between the state m s — +1 and m s = — 1 in the 
large detuning regime, (c) The effective coupling provides a 
dressed state basis which then permits energy conserving flip- 
flops between the dressed states and external resonant nuclear 
spins. 


II. MAGNETIC MANIPULATION AND 
POLARIZATION OF A SINGLE NUCLEAR SPIN 
NEARBY A RANDOMLY ORIENTED NV SPIN 

A. Difficulties from the random orientations 

The nanodiamond ensembles that we are interested in 
are realised as powder or solutions. In both cases the 
spatial orientations of the nanodiamonds and therefore 
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of the NV centers hosted in them is random and uni¬ 
formly distributed across the full solid angle (see Fig. la 
for an illustration). In the laboratory frame this leads 
to the lack of a common quantization axis for the NV 
centers. We begin by describing in more detail the two 
principal challenges that are being imposed by these ran¬ 
dom spatial orientations: (1) The direction of the natural 
quantization axis associated with the crystal-field energy 
splitting D is not controllable, resulting in a significant 
variance of the energy levels in the presence of an ex¬ 
ternal magnetic field and (2) optical pumping of the NV 
center electron spins in an external magnetic field will 
not initialize all NV center electron spins to the same 
state. 


1. Zero-field and external magnetic field distribution 

The negatively charged NV center, in the following de¬ 
noted for brevity as the NV center, realises in its ground 
state an electronic spin triplet (5 = 1 ) which exhibits 
a zero-field splitting of D = (27r)2.87 GHz which sepa¬ 
rates the \m s = 0 ) state energetically from the degenerate 
\m s = ±) manifold. The application of an external mag¬ 
netic field, in the following assumed to take the value 
B = 0.36 T, lifts this remaining degeneracy such that 
the state m s = — 1 is shifted below the m s = 0 state (see 
Fig. [TJd). The zero-field splitting and the Zeemann effect 
due to an external magnetic field of the NV center are 
described by 


H NV = SD§+'y e BS. (1) 

Herein D denotes the orientation dependent zero-field 
splitting tensor, y e = (27r)28.7GHz/T the gyromagnetic 
ratio, S the electron spin -1 vector operator [33] and B 
the magnetic field vector. In the principal axis system 
defined by the NV symmetry axes, the zero-field splitting 
tensor D is diagonal 

D = diag(-^D +E,~D - E,^D). ( 2 ) 

Here E denotes the strain dependent contribution. It 
is worth emphasizing that in the scenario that we are 
considering here, the orientation of the symmetry axis 
of the NV center relative to the external magnetic field 
is uniformly distributed over the unit sphere. For the 
following it will be convenient to conduct the discussion 
in the laboratory frame whose z-axis we define to take the 
direction of the externally applied strong magnetic field, 
7 e B D. In this frame, the zero-field splitting tensor D 
will have off-diagonal elements. As 7 e B ^ D these off- 
diagonal elements are rapidly rotating so that their main 
effect will be energy shifts of the diagonal element of the 
Hamiltonian. The Hamiltonian can then be written as 
(see the Appendix for details of the derivation) 

H'J ff = ('leB + 6(9))S Z + D(8)S 2 Z . (3) 
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FIG. 2. a, Zero-field distribution D(0) of the NV spins in nan¬ 
odiamonds with D = (27r)2.87 GHz and E = (27t)20 MHz. b, 
The second order corrections which induce a energy distribu¬ 
tion 5(0). 


where 


D{0) = T>(1 + 3cos(2$)) + 3E(1 — cos(2$)) ^ (4) 


m = 


7e£|G! 


+ 


|G 2 


{leB) 2 — [D(0)] 2 2 7e H’ 


6 is the angle between the magnetic field direction and the 
NV-axis and the Gi are given in the Appendix. Clearly, 
the random orientations of the NV centers cause a vari¬ 
ation of the zero-field splitting D(0) across the entire in¬ 
terval [— (27t)1.43GHz, (27t)2.87GHz] and 5(0) across the 
interval [0MHz, (27t)140MHz] as shown in Fig. [2j 

If we were to follow the scheme of Ref. [24| , i.e. without 
the combination of off-resonant drive and ISE technique 
that we will present in this work, the uncertain detuning 
of the MW frequency from the electronic resonance can 
reach the order of GHz, which in turn would prevent 
effective polarization transfer. 


2. Optical initialization of NV spins at the high magnetic 
field limit 

A second important difference between a randomly ori¬ 
ented nanodiamond ensemble and a bulk diamond con¬ 
cerns the optical polarization of electron spins of the NV 
center. For bulk diamonds, the magnetic field can be 
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aligned with the principal axis of the NV center and the 
electronic spin of the NV center can be optically polar¬ 
ized to the state \m s = 0) by illumination with a 532 
nm green laser. However, for an ensemble of randomly 
oriented nanodiamonds, even in the limit of a strong mag¬ 
netic field, the NV centers will be optically pumped to the 
state \m s = 0}e that is defined by the relative orientation 
of the NV center with respect to the externally applied 
magnetic field which defines the laboratory frame. 

As discussed in Appendix, these two coordinate sys¬ 
tems can be transformed into each other and, employing 
of S Ze = cos 0S Z — sin #(cos 05^ — sin (j)S y ), we can ex¬ 
press the eigenstate \m s = 0 )# in terms of the eigenstates 
|0), | =b 1} of eq. (|3]), i.e. in the lab frame, as 

|0 > 9 = cos(91 0 ) + ^(e*| + 1 ) - e-*| - 1 ». (5) 

If 0 is large, the eigenstate |0) of the NV center in the 
laboratory frame as given by Hamiltonian eq. © differs 
significantly from the zero-field eigenstate | 0)<9 of the NV 
center. Hence optical initialization of randomly oriented 
NV centers lead to very different states depending on the 
orientation of the NV center. 

Notice though that for moderate misalignment be¬ 
tween the NV center and external magnetic field, the 
initialisation of the NV center is well approximated by 
|0). Indeed, for 0 < 10° degrees, |(0|0 )^| 2 > 0.97 and 
for 0 < 20 ° we find |( 0 | 0)#| 2 > 0.88 which implies sig¬ 
nificant polarization along the quantization axis defined 
by the external magnetic field. Therefore, for all the ori¬ 
entations of the NV centers that fall into two spherical 
sectors whose cone angle is 2 6 , i.e., the two blue spheri¬ 
cal sectors in Fig. [IJi, significant optical polarization can 
be achieved. From now on, we define a deviation by 0 
to imply that the symmetry axis of the NV spins falls 
within these two spherical cones. 


B. The main ideas 

The main idea that we will develop in this section con¬ 
sists of three main ingredients, namely (i) the genera¬ 
tion of an energy gap between electronic states that is 
relatively robust with respect to variations in the rela¬ 
tive orientation between the NV center and an external 
magnetic field and at the same time (ii) a strong cou¬ 
pling between these electronic states and nuclear spins. 
Finally, (iii) the integrated solid effect is employed to 
achieve additional robustness against imperfections. All 
this is achieved by the use of dressed states in a ladder 
type configuration in which the single quantum transi¬ 
tions are far detuned and the double quantum transition 
is nearly resonant (see Fig. EJ). 

In such a setting, that is for a detuning on the single 
quantum transition that is A(0), the effective Rabi fre¬ 
quency and thus the energy splitting on the double quan¬ 
tum transition is proportional to Vt 2 / A (0). In leading 


order in 0 , this suffers a variation that is of the order of 
Q 2 /A 2 This contrasts with the case of detuned driv¬ 
ing one single quantum transition only. In that c ase the 
effective energy gap is proportional to VA 2 + Q 2 which, 
if A >> £2, suffers a variation that is of order |A Q m There¬ 
fore the variation in the former is suppressed by a factor 
of Q 2 / A 2 compared to the latter. This makes it much 
easier to achieve and maintain a Hartmann-Hahn reso¬ 
nance. 

The second crucial ingredient is related to the fact 
that the coupling between the electron spin and the nu¬ 
clear spins is mediated by terms of the form S z ® I z . 
In a far-detuned single quantum transition the eigen¬ 
states involved in the Hartmann-Hahn resonance are 
| m s = 0) + £|m s = + 1 ) and \m s = + 1 ) - £| m s = 0). 
As these states are close to being eigenstates of S z they 
hardly couple and in consequence the electron-nuclear 
coupling is suppressed by a factor proportional to £ 
when compared to the coupling between states equally 
weighted superpositions \m s = 0) zb \m s = +1) which 
one would have obtained for resonant driving. Crucially, 
thanks to the small detuning on the double quantum 
transition we obtain \m s = — 1 ) zb \m s = + 1 ) which lead 
to a strong coupling between electron and nuclear spins if 
a Hartmann-Hahn resonance is realised (see r.h.s of Fig. 
3). The latter is relatively robust thanks to the strong 
detuning on the single quantum transitions. 



FIG. 3. The level structure of the three level Hartmann-Hann. 
A gap is created between the | D) = -^= (| + 1) + | — 1)) and 

the | B) = -j= (| + 1) — | — 1)) states. This gap is robust to 
changes in the one photon detuning (A) and is also robust 
to small changes in the two photon detuning((5). Due to 
the fact that the coupling to the nuclei is via the dressed 
states(|D), |B)) and not the bare states the coupling is not 
decreased due to the off resonant drivings. 

So far we have ignored additional energy shifts on the 
double quantum transition 5(0) which will lead to a loss 
of the Hartmann-Hahn resonance. In order to confer ad¬ 
ditional robustness to our scheme, rather than aiming to 
maintain a fixed Hartmann-Hahn resonance, we consider 
an adiabatic sweep of the external magnetic field on the 
double quantum transition where start and endpoint of 
the sweep are chosen such that intermittent Hartmann- 
Hahn resonances are guaranteed. This technique, the 
integrated solid effect, allows polarization transfer from 
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a large fraction of the NV electron spins to their sur¬ 
rounding nuclear spins. 

Therefore, the resonant double quantum transition 
that is mediated by far-detuned single quantum tran¬ 
sitions combines robustness with strong electron-nuclear 
coupling. It is this idea and its combination with the inte¬ 
grated solid effect that we are going to develop and anal¬ 
yse in detail for realistic experimental parameters and 
various imperfections in this section. 


1. Double quantum transition 


For illustration of the basic idea, let us first consider a 
system that is composed of an NV center and a single 13 C 
nuclear spin in a nanodiamond of random orientation. In 
the laboratory frame the applied magnetic field defines 
the z-axis, the NV is placed at the origin of the coordinate 
system and the 13 C nuclear spin, placed near the NV 
center, is situated at position r. A microwave (MW) 
field of frequency ujm is applied as off-resonant drive of 
the spin transitions | — 1 ) | 0 ) and | 0 ) | + 1 ) while 

approximately satisfying a double resonance condition, 
that is 2ujm — E\+)—E\_y We chose circular polarization 
but linear polarization would suffice. The Hamiltonian of 
the whole system is then 


H = cos uo M t + S y sin u M t) + ( 7 e E + S(0))S z 

+D(0)S 2 +g\si-3(S-e r )(I-e r )1 + ln BI z , (6) 


in which g = with r = |r| denoting the distance 

from the NV spin to the nuclear spin and e r = f/r. I is 
spin- \ vector operator of the nuclear spin with the gyro- 
magnetic ratio y n and Qm = a/ 2 ^ is the Rabi frequency 
of the driving field. Assuming a point-dipole interaction 
and neglecting the contact term we obtain 


H e ~ &m{S x cos UMt + S y sin wm^) 

+{ le B + 5 (0))S Z + D(d)S 2 z + ln BI z 
-gS z [3 e z r (e x r I x + el>I y ) + (3(e ):) 2 - l)I z 


(7) 


In the interaction picture with respect to Ho = ujmS z , 
the Hamiltonian is given by 

H' = Q(\ — 1)(0| + 10) (+11 + h.c.) 

+ ln BI z + D(9)S 2 + AS Z + S z ■ A ■ I, ( 8 ) 

in which A = g\J 1 + 3(e*) 2 h = Ah is the hyperfine inter- 
action tensor for the nuclear spin, with h determined by 

e r : h x = 3e^eV v /l + 3( e^) 2 , h v = SeJ'e'/y'l + 3(e):) 2 , 
and h z = (3(e ^) 2 — 1 )/ \/l + 3(e,;') 2 . The detuning A is 
given by 


A — 7 e B + 5(0) — cjm, (9) 

which is a function of the NV orientation due to the sec¬ 
ond order energy correction 5(0). 


Since the nuclear coordinate system is arbitrary, we 
can redefine it by following the method in Ref. [24[. The 
z-axis is chosen to coincide with the external magnetic 
field, and we define the y-axis to be perpendicular to 
both magnetic field direction and h. In the new nuclear 
coordinate system, the Hamiltonian is given by 

H" = fi(| - 1)(0| + |0)(+1| + h.c.) + D(9)S 2 + AS Z 
AlnBIz’ + S z ■ ( a x 'l x i + a z 'l z /) (10) 

where I z = I z r and a Z f and a x > are the elements of the 
secular and pseudosecular hyperfine interactions, respec¬ 
tively. 

In what follows we assume that (i) \D(0)\ A which 
implies that the microwave driving field tends to be far 
detuned from the single quantum transitions | — 1 ) | 0 ) 

and | 0 ) | + 1 ) and (ii) \D(6)\ ft such that the 

Rabi frequency is significantly smaller than the detuning 
(see Fig. lib). This will allow us to remove level |0) 
adiabatically from the dynamics to achieve a simplified 
effective Hamiltonian eq. (HU). 

To this end, we consider the dominant parts of eq. 
m by neglecting the nuclear spins as well as A S z . In 
the basis {| + 1 ), | 0 ), | — 1 )} we therefore consider 

( d{ 0) n o \ 

h nv = O 0 0 . (11) 

\ o n d{9) ) 

The eigenstates and eigenenergies of are 

, 1 (| + 1) + |-1)-X ± |0)), (12) 

V 2 +^± 

= \ [D{9) ± y/ 8 fi 2 + D 2 (0)] (13) 

and 

|A) = +(| + 1 )-|- 1 )), (14) 

= D[9) (15) 

with 


For \D(6)\ > 9 we have two cases: D(0) > 0 leads to 
X-(D(0),Q) - 0 and X+(D(0),Q) > 1, while D{0) < 0 
results in X+(D(6),Q) ~ 0 and X-(D(6),Q) 1. 

In both cases, two of the eigenstates are approximately 
given by |±) = -^=(| — 1 ) ± | + 1 )) and form an effective 
two-level system while the state | 0 ) does not participate 
in the dynamics because for \D(0)\ O it is far detuned. 
We will now focus attention on these two states |±) and 
write the Hamiltonian eq. m in the subspace spanned 
by (|+), |-)} to find 

— di Cl e ffU z T 2A (j x T , y n BI z / -t- ‘2(j x )a x ' I x > -t- ci Z 'I z /) 

( 16 ) 
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in which H+ and H_ corresponding to the cases D(9) > 0 
and D(0) < 0, respectively. Here 

Veff = \[-\D{6)\ + y/8WTmO)} 

and a, = ±(|+)<+| - |-)(-|), = ±(|+)<-| + |-)(+|), 

Iz> = |(l txt I - I 1)4 I) and I x , = I(| t)4 I + I m !)• 
Furthermore | and | f) denote the ground and excited 
states of the nuclear spin. 

For the final step, we now restrict attention to the 
case where 9 G [0°,20°] in which case D(9) > 0 and 
we consider Hamiltonian H+. In matrix notation, the 
electronic part of the Hamiltonian eq. m then takes 
the form 



The eigenstates and eigenenergies are 

lx-) =cos||-) -sin||+) 

|X+) = cos^|+) +sin||-), 

CVeff = ± ^ A 2 + —(18) 

in which arctan£ = qA^/2 •> ^ and ^eff are dependent 
on the angle 9 between NV orientations and applied mag¬ 
netic field. Defining Pauli operators 07 , 07 , erg in the 
eigenbasis of H& we find 

H+ = 2 LdeffCTz + 7 nBI z ' (19) 

+ 2(07 sin p + 07 cos cp) • (« a x >I x > + a Z 'I z i), 

in which the angle <p is determined by 



The Hamiltonian H+ can be simplified further under 
two main assumptions, namely (i) a x > <C 7 n F>, that is for 
weakly coupled nuclear spins and (ii) | 2^ — uj e f /1 « a x > 
that is we satisfy a Hartmann-Hahn condition. Then the 
polarization transfer dynamics is described by 

Htrans = 2(J e ff^z + lnBI z ' + 2 a z > COS p(J z I Z ’ 

a T > sing? A , 

H- 2 -(lx+jl)(X-?t I + (20) 

where 2a z > cos pcr z I z ' does not affect the flip-flop interac¬ 
tion. 

Let us briefly highlight the key difference of this off- 
resonant driving in comparison to previously introduced 
polarization schemes, namely its relatively high robust¬ 
ness with respect to the angle 9 between the external 
magnetic field and the axis of the NV center. If we were 
to follow the scheme of Ref. [24[, even for 0 = 1°, the 


detuning of the MW frequency from the NV electronic 
resonance exceeds (27r)l MHz, which prevents efficient 
polarization transfer (for 6 = 20° we find (27r)500 MHz 
detuning). In our scheme, for as long as D(6) D, 
the energy difference between the states |—) and |+) 
scales as 2 Vt 2 /D(0) and so does the effective Rabi fre¬ 
quency Cl e ff- For deviation 6 G [0°,20°], if D = (27r)65 
MHz, the effective Rabi frequency is roughly within 
[(2tt) 3MHz, (2 tt) 3.6MHz], which narrows the range of 
the detuning of the resonant frequency. 


2. Integrated solid effect 

In the discussion so far, we have assumed that the 
electron spin of the NV center is held continuously at 
a Hartmann-Hahn resonance with a specific target nu¬ 
clear spin to achieve polarization transfer. In practice, 
however, this will only capture a small fraction of the 
NV centers and nuclear spins as the resonance condition 
will depend on A **= 7 e B + 5(0) — ujm which in turn is a 
function of the angle 9 via 5(9) (Fig. [ 2 ) 3 ). In the range 
9 G [0°,20°], we have 5(9) G [0, (2tt)45MHz], as shown 
in Fig. [2b which can lead to a violation of the Hartmann- 
Hahn resonance for a large fraction of the NV centers. In 
order to achieve polarization transfer for a larger fraction 
of the NV centers, we will use our ability to control the 
strength of the externally applied magnetic field or the 
frequency of the applied microwave field to implement a 
sweep of A. 

In order for us to understand qualitatively the condi¬ 
tions that such a sweep has to satisfy, we briefly examine 
the energy level diagram of a system consisting of one 
electron and one nuclear spin as described by eq. ([ 20 ]). 
This reveals that there are two resonance points A\ and 
A 2 for A = ±A hh owing to the symmetry of the Hamil¬ 
tonian. In a fully adiabatic sweep, polarization would 
exchange twice and result in a vanishing net polarization 
transfer. Hence the sweep rate needs to be chosen such 
that around the A{ the system becomes non-adiabatic for 
the two branches including |x+, |) and \x~A) as this will 
lead to an approximately 50% probability for polariza¬ 
tion transfer. At the same time the sweep has to remain 
adiabatic with respect to the energy gap for transitions 
between the {|x+4), |X->t)} and the {|x+,t). IX-4)} 
manifold as transitions between these two manifolds, in¬ 
duced by non-adiabaticity or in fact dephasing events, 
would lead to depolarization. These requirements set 
some limitations on the sweep rate that are, fortunately, 
not too stringent. 

The quasi-adiabatic sweep described above realises an 
instance of the so-called integrated solid effect (ISE) 
and can achieve polarization transfer between a large 
fraction of the NV spins and nuclear spins. Therefore, 
magnetic control of the nuclear spin, and efficient po¬ 
larization exchange, become possible even for relatively 
large angles between NV-center axis and external mag¬ 
netic field. In the following we will examine this idea in 
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FIG. 4. Eigenenergies of eq. (l20l) for B m 0.36 T in which 
| a * /S 2 in( ^ | = (2 tt)0.5 MHz, Q eff = (2tt)2.2 MHz (the dot- 
dashed red line), H e // = (27r)3 MHz (the dotted green line), 
Q e ff — (27r)3.6 MHz (the dashed blue line). The level dia¬ 
gram shows the progression of the eigenvalues of the Hamilto¬ 
nian during the ISE sweep. Energies of the states |x+,t) and 
|X-,f) (solid blue line) for Cl e ff — (27r)2.2 MHz are included. 
Transitions from states represented by solid blue lines to state 
represented by dashed lines may lead to depolarization of the 
nuclear bath. 


and indeed during the entire sweep there will be no tran¬ 
sitions of the type |x-,i) to |x+,i) which, as is easily 
seen by examination of Fig. 01 would have a depolarizing 
effect in the subsequent sweep. 

To gain a feeling for the upper limit on the sweep rates 
that this implies for typical experimental parameters sup¬ 
pose as in the previous section that Q = (27r)65 MHz 
so that in the range 6 G [0°,20°] we find fl e ff m the 
range [(2tt)3MHz, (2tt)3.6MHz\. This in turn implies 
that v < (27r)10 MHz//is guarantees that the adiabatic 
condition is satisfied. As shown in Fig. [2t> in this range 
we have 5(0) G [0, (27r)45 MHz\. Suppose a sweep rate 
v = (27t) 6 MHz//is, then no more than 10 /is are required 
to cover all NV centers within 0 G [0°,20°]. Experimen¬ 
tally such sweep rates are easily obtained using arbitrary 
wave form generators. 

Let us now proceed to discuss case (ii), namely the 
sweep around the Hartmann-Hahn resonances which will 
provide lower bounds in the sweep rate. Near the res¬ 
onant points, the effective Hamiltonian eq. ([20]) in the 
subspace spanned by {|x+, |), |x-? t)} is well approxi¬ 
mated by the matrix form 


more detail. 

Let us now examine two regimes of the quasi-adiabatic 
transfer described above, namely (i) those parts of the 
sweep parameters that are far from the Hartmann-Hahn 
resonance points A\ and A 2 which will lead to an upper 
bound on the sweep rate and (ii) the behaviour around 
the Hartmann-Hahn resonances which will provide lower 
bounds on the sweep rate. 

For regime (i), as the coupling between electron and 
nuclear spin is negligible, it suffices to examine the elec¬ 
tronic part of eq. cni, that is 


it’ _ I ~2~ *(<> 

4 “ l A«) - V 


( 21 ) 


in which A (t) = 7 e B — 5(6) — ujm • For simplicity we 
assume the detuning to vary at a constant rate v > 0 , i.e. 
A (t) = A (ti) + vt , where A (ti) is the initial detuning at 
the start of sweep. According to the quantum adiabatic 
condition we require, at all times, that 


(El m (t)\Ei n (t)) 
\Elm{t) ~ E ln (t)\ 


m V n, 


( 22 ) 


where E\ m (t) and E\ n (t) are the instantaneous eigenval¬ 
ues, \Eim(t)) and \Ei m (t)) are the eigenstates of effec¬ 
tive Hamiltonian respectively. Then the condition 
for the sweep to be adiabatic is given by 


(-l) l v(t-ti) cos ipi ax ' s 2 in - 

— 2 ° - (~l) l+1 v(t -ti) COS (pi 

(24) 

Here Z = 1 , 2 corresponds to points A\ and A 2 , and near 
these two points we make use of the expansion 

^eff(t) = Veff(tl)+ dU)e g^ t '> | = (-l) l v(t-ti) COS (fit. 

According to Landau-Zener (LZ) theory [35|, transitions 
are possible between two approaching levels as a control 
parameter is swept across the point of minimum energy 
splitting. The asymptotic probability of a LZ-tunneling 
transition is given by 



Plz = e 




(25) 


where, by virtue of the Hartmann-Hahn condition 7 n B = 

y^4A 2 + ^eff resonan t pints A\ and A 2 , we find 


M = 


a x f sin (pi \ 2 

2 / 


Q 2 


e// a ^' 


2M c osw 8|«|(7 n B)y/('y n B) 2 -Q 2 eJf ' 

This enters the unitary transformation 


N L = 


VI - PlzPV* 

VPlz 


~VPlz 

VI - PLze-^s 


(26) 


If » I- (23) 

Adiabatic evolution implies that the initial eigenstate of 
this Hamiltonian H/± will remain close to the instanta¬ 
neous eigenstate at any time during the sweep. In partic¬ 
ular, during approach to the Hartmann-Hahn resonances 


that is taking place when traversing the avoided crossing. 
Here (p s = cp s — 7 r /2 is due to the Stokes phase ip s . 

The full sweep takes the system from A <C 0 to A >> 0 
across both A\ and A 2 . At each resonance point po¬ 
larization transfer is possible between the electron and 
the nuclear spins J36|. Consider a sweep starting at 
A( 0 ) <C —&eff fc> r which the electron spin is initialized 
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in state |y_) and the adiabatic condition eq. (|23|) is satis¬ 
fied. If the nuclear spin is initially found in state | J,), the 
electron-nuclear system will start and finish the sweep in 
state |y_,^), i.e. the nuclear polarization is unchanged. 
If on the other hand the nuclear spin is initially found in 
state | t), the evolution of the system will start in |x_, t)- 
In this case the probability for experiencing a polariza¬ 
tion transfer at the end of the full sweep, that is finding 
the system in state |y + ,|), is given by 

P = Pmax Sin 

= 4Plz(1 - Plz) sin 2 [4>s't]. (27) 

Here &st is determined by the initial phase of the system 
and the phases acquired during the adiabatic evolution 
and the non-adiabatic transitions. We do not discuss 
these dynamical phase shifts in detail as they depend on 
details of the evolution and orientations that vary from 
NV center to NV center and it is therefore reasonable to 
assume that &st is random [33, 38]. Hence, when aver¬ 
aged over many sweeps, the average polarization transfer 
probability is given by 


P = 2Plz(1 — Plz), 

as shown in Fig. [5j The value of Plz and hence the 
maximum P = 2Plz{^ ~ Plz ) is determined by the hy- 
perfine coupling strength a x /, the rate of the ISE sweep 
v, and to a lesser extent on the strength of the applied 
microwave field. 

It is straightforward to see that Pmax and hence P 
vanish for both v 0 and v —>> 0, corresponding to the 
diabatic and adiabatic sweep. The maximum value is 
achieved for Plz = 1/2 and the parameters for achieving 
this maximum are quite flexible, see Fig. [6] Pmax re¬ 
mains high across a wide range of sweep rates as shown in 
Fig. [6^ and b. Note furthermore, that for a given sweep 


Energy 



FIG. 5. Double-passage transition. A schematics of a pair 
of the adiabatic energy levels as in Fig. [H At each crossing 
there is a probability of Plz for a Landau-Zener transitions. 
The lines with one (two) arrows show the two paths where the 
transition to the upper level happens during passage of the 
first (second) resonance point. The average probability for a 
polarization transfer at the end of the completed passage is 
given by 2Plz(1 — Plz)- 


rate we obtain efficient polarization transfer to the nu¬ 
clear spins in a broad range of different microwave Rabi 
frequencies, see Fig. [6l For example, for ft e ff = (27t)3.6 
MHz, a speed v = (27t)6 MHz//is can obtain Pmax > 0-2 
in the coupling range a x > G [(27t)0.18, (27t)0.9]MHz. Fur¬ 
thermore, a slower sweep gives a better polarization 
transfer for a smaller effective Rabi frequency and a 
weaker coupling strength between the electron spin and 
the nuclear spin. In order to polarize the nuclear spins 
with larger distance from the NV spin, one can slow down 
the ISE sweep, or increase the Rabi frequency of the mi¬ 
crowave field, as shown in Fig. [6t. 


C. Preparation of the initial state 


So far we have assumed the preparation of a specific 
initial state in the polarization sequence to be achieved 
with perfect fidelity. However, due to the random ori¬ 
entations of the NV centers, optical pumping leads to 
the initialization of the electron spin of the NV center 
in a wide range of different states. Furthermore, due 
to the broad distribution of zero-field splittings of the 
NV center in randomly oriented nanodiamonds, any ap¬ 
plied microwave field will experience an uncertain and 
potentially large detuning from the electronic resonance. 
This presents challenges to the preparation of the elec¬ 
tron spin of the NV center in the initial state that is 
required for our polarization scheme. We address these 
challenges with an adiabatic sweep of the frequency of 
the microwave field which we will show to deliver robust 
and rapid state preparation. 

Let us consider as an example the initialization of the 
NV center in the state | — 1). For moderate 6 optical 
pumping will lead to a state that possesses a significant 
overlap with |0). In order to map this state to the target 
| — 1), we make use of the fact that in the relatively 
high magnetic fields that we are considering (e.g. B= 
0.36T) the energy gap on the | — 1) |0) transition 

is very different from the energy gap of the |0) O | + 1) 
transition. Hence it is possible to use a circular polarized 
microwave driving field 

Hdr = ^m-{S x cos uj-t + S y sincj_t), (28) 

in which Qm- = \/2f2_ is the Rabi frequencies of the 
MW field, with frequency to drive dominantly the 
| — 1) |0) transition. The resulting effective two-level 

system consisting of the states |0) and | — 1) can be de¬ 
scribed in the matrix representation by 


( A M w(t)/ 2 \ 

\ £2- — AmwW/2 ) ’ 


(29) 


in which A Mw(t) = + 5(0) + D(6) — u-(t). Now 

we assume that uj-{t) experiences a constant rate sweep 
such that A Mw(t) = A Mw(ti) + A_(£ — +). We set 
U = 0 for simplicity. The condition for the sweep to be 
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FIG. 6. The maximal polarization transfer P macc as a function of the effective coupling Q e //- a, Q e ff — (27r)3.6 MHz b, 
Q e ff = (27r)3 MHz; c, Q e ff = (27r)2.3 MHz. 


adiabatic is then 


D. Polarizing cycles and effective ranges 


8Q 2 


> i. 


(30) 


Adiabatic evolution implies that the initial eigenstate of 
this Hamiltonian i^o,-i will remain close to the instan¬ 
taneous eigenstate at any time during the sweep. The 
eigenstates and eigenenergies are 

l^ ± W) = Ci T (i)|-l)±Ci ± (i)|0), 

E ± (t) = (31) 


with E+ E+(t) M ^ ^ a ^ the s ^ ar t °f the sweep 

Avpm( 0) < 0 and \Awm(0)\ 2f2_ and the electron 

spin is prepared in the state |0), the state will evolve 
along the path \E+(t)} to end in the state | — 1). 

As an example, consider a Rabi frequency = (27 t) 20 
MHz and 6 G [0°,20°]. Then a sweep from A wm(U) > 
(27t)160 MHz to A< — (27 t)160 MHz includes all 
NV centers in the sweep. For | — u)\(ti)\ = uj\tf ~ 
(27t)870 MHz we require tf [uj\{tf) — uj\(ti)\/8Pt 2 _ = 
0.04 //s to achieve adiabaticity. Hence very high fidelity 
preparation of the state | — 1) can be achieved within 0.4 
fi s which is faster than the optical polarization cycle. 

So far we have assumed that the optical initial state 
|0) is well prepared. However, for the relatively high 
magnetic fields that we are considering, the optical initial 
state is given by |0)^ = cos 0|O) + + 1) — e~ lcf) \ — 

1))], (see eq. (|5j)). For a perfect adiabatic sweep we 
obtain the average polarization in the range 0 G [0°, 20°] 
as 



Here S is the surface area of the two sectors within 6 
deviation and dS is the corresponding area element. The 
reachable polarization for nearby nuclear spins and the 
average polarizations are P^ 1 = 0.91 for 6 G [0°, 20°]. 


We can now summarize the polarizing cycles as follows 
(see also the black box of the flow diagram Fig. ( 3 ): (i) 
Optical pumping initializes the NV spin and (ii) a sub¬ 
sequent adiabatic transfer brings the electron spin to the 
| — 1) state, (iii) By using off-resonant driving and an 
quasi-adiabatic sweep across two Hartmann-Hahn reso¬ 
nance points (ISE technique) polarization is transferred 
from the NV spin to the nuclear spin. The steps (i)-(iii) 
are repeated as required. This method has two primary 



Adiabatic 

rotation 



Polarizing cycles 


FIG. 7. A simple flow diagram of the entire polarization 
scheme: For a given direction of magnetic field, the large 
box includes numerous polarization cycles (each polarization 
cycle includes optical initialization of the NV spin, adiabatic 
preparation of the state | — 1) and then polarization transfer 
by matching a H-H resonance condition by using off-resonant 
microwave driving and the ISE technique). After the adi¬ 
abatic rotation of the magnetic field (The thick black arrow 
represent the magnetic field orientation which is rotated to an¬ 
other direction denoted by the dashed blue arrow) we repeat 
our polarizing cycles to achieve the polarization of additional 
nuclear spins. 
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FIG. 8. Exact polarization dynamics with B — 0.36 T, v = (27r)6 MHz / ps and At = 10 /xs: a, The coupling strength is given 
as a x f = (27r)0.6 MHz, U e // = (27r)3 MHz (blue circles) and Q e ff — (27r)3.6 MHz (red triangles), b, We have U e // = (27r)3.5 
MHz, a x f — (27r)0.7 MHz (red triangles), a x > — (27r)0.5 MHz (blue circles) and a x > — (27r)0.3 MHz (green cubes), c, Exact 
polarization dynamics with v = (27r)0.8 MHz//xs, At = 50 /xs, and Q e ff = (27r)2.3 MHz, a x r = (27r)0.3 MHz (red triangles) 
and a x f = (27r)0.1 MHz (blue circles). 


advantages: (1) the impact of the random orientations 
is minimized, and NV spins in a wide range of orienta¬ 
tions are still accessible for polarization. (2) Using this 
method, the coupling strength between the NV and nu¬ 
clear spins is doubled because it involves the double quan¬ 
tum transitions between the states | — 1) and | + 1). 

In Fig. [8] we simulate many cycles of this polarization 
protocol for a system of one NV center and a single nu¬ 
clear spin under the assumption that the NV center is ini¬ 
tialized to the state | — 1) as an example. The initial den¬ 
sity matrix of the unpolarized nuclear spin is po = 1 / 2 , 
with I being the unit matrix. For the state pk of the 
nuclear spin after the k-th iteration of the polarization 
protocol we have the following relation 

Pk +1 = tr e [U t {p k <g> \x-)(x-\) u t] (33) 

where U t = e ~ f lHtransdt j s time evolution operator, 
Htrans is the Hamiltonian eq. flU, we denote by At the 
sweep time and by tr e the trace over the electron spin. 
The polarizing P is defined as P = (I z f )/(/*')o with (I z r) 
denoting the expectation value of the nuclear spin and 
(I z ')o the expectation value of the completely polarized 
state. 

Finally, as we have discussed in the previous sec¬ 
tions, the ISE requires an adiabatic slow passage. For 
9 G [0°,20°] we find 5(6) G [0, (2ir)4:bMHz] and in order 
to include the entire range of detunings in the sweep we 
can estimate the required polarizing time. At a sweep 
rate of v = (27r)6 MHz / ps we find At ~ 10 /xs. In Fig. 
[8fei, assuming the perfect initialization of the NV spin in 
the state | — 1), we take one nuclear spin with the typical 
coupling a x > = (27r)0.6 MHz as an example to demon¬ 
strate the polarization transfer (the nuclear spin is about 
0.5 nm from the NV spin and a z r = (27r)0.64 MHz). Here 
the polarization transfer after the first cycle, Pi, does not 
quite achieve the maximum possible value in Fig. [6] due 
to the existence of the phase shift &su more precisely 
Pi = 0.58 < Pmax = 0.7 when Q e ff = (27r)3 MHz, see 
both Fig. [6 Jd and [H^t. Simulations for different hyperfine 
couplings are shown in Fig. [HJd. Again the sweep rate 


v = (27r)6 MHz / ps and polarization is transferred to 
the nuclear spins with different coupling to the NV spin, 
which agree with Fig. [6j Additionally, slowing down 
the ISE sweep, as shown in Fig. [8t for a sweep rate 
v = (27r)0.8 MHz / ps we can polarize more distant nu¬ 
clear spins with a coupling strength a x > = (27t)0.1 MHz, 
which corresponds to a distance of about 1 nm from the 
NV center. This coupling is approaching the limitation 
of our polarization scheme, since the rotating-frame spin- 
lattice relaxation time of NV spin in nanodiamonds is 
estimated to be about 100 /xs which places a lower limit 
of the duration of the full sweep [39]. 


III. THE BENEFIT OF ADIABATIC 
ROTATIONS 

The NV polarization in our scheme is limited by the 
adiabatic sweep range and efficiency of the initialization 
of NV spins, enabling the polarization up to a maxi¬ 
mum angle between NV center and external magnetic 
field. As discussed above, our scheme certainly performs 
well for NV centers whose orientation forms an angle of 
less than 20 degrees from the externally applied magnetic 
field (and may potentially work for somewhat larger an¬ 
gles). It is however desirable to extend the polarization 
to a larger fraction of the nanodiamond ensemble, ideally 
covering the entire solid angle. One method for achiev¬ 
ing this applies an adiabatic change of the orientation 
of the magnetic field relative to the NV center, either 
by rotation of the magnetic field or by rotation of the 
nanodiamond. For example, following the magnetic field 
rotation to the Z'-direction, another set of NV centers 
with orientations in a solid angle (shown in red in Fig. [9]) 
around the new magnetic field direction will participate 
in the polarization dynamics while the NV centers in the 
original orientations (the blue cone in Fig. [9j) will be in¬ 
active. Hence previously polarized nuclear spins will not 
affected by the polarization sequences, but these nuclear 
spins may use this ’’idle” time to spread their polariza- 
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Step2 


Step3 



Polarizing cycles 



FIG. 9. a, Diagram of adiabatic rotating direction of the magnetic field. Blue spherical sectors present 0 degrees deviation 
tolerance for Z-direction of the magnetic field. Magnetic field is rotated by angle $ to get red spherical sectors with the half¬ 
angle of the cone angle O' to be involved, b, Flow diagrams of polarizing steps of our scheme related to the adiabatic rotations 
of the magnetic field. The small blue and red arrows are the nuclear spins nearby the different NV spins in different spherical 
sectors. Four steps are given for our polarizing scheme. We use 2D plotting for simplicity, the sphere is presented by a circle 
and the areas formed by the dashed lines denote the corresponding spherical sectors. It is shown by the arrows that all the 
nuclear spins follows the rotated magnetic fields. 


tion by spin diffusion across larger volumes until further 
rotations subject them to another polarization cycle. 

More specifically, let us consider the rotation of the 
magnetic field as an example. The Hamiltonian that is 
describing the rotation of the magnetic field from the Z- 
direction to the Z'-direction is given by 

Hb = 7atHz[cos(^£)/z + sin (M)Ix] = 

where $(t) = fit is the angle between the rotated mag¬ 
netic field and its initial direction and 7 nBz = (27t) 4 
MHz. 

Similar to the discussion in the initialization scheme, 
the well-known adiabatic condition requires $ 2 yA tBz- 
For a rotation by $ = 180° this leads to the condi¬ 
tion t r 0.078 /is for the minimal rotation time t r . 
This implies that t r = 0.78 /is will ensure adiabaticity 
and minimal perturbation of nuclear spin polarization 
for a rotation to an arbitrary orientation. This rotation 
time is much shorter than even a single polarization cy¬ 
cle. Hence, after running a few polarization cycles that 
achieve a high net polarization of nuclear spins for NV 
centers oriented along the Z direction, we can adiabati- 
cally rotate the magnetic field to another direction. The 
polarization of the nuclear spins will then be achieved for 
this new direction, as shown in Fig. [9] Note though that 
the rapid rotation of a strong magnetic field is challeng¬ 
ing. Hence one needs to resort either to a mechanical 
rotation of the sample as a whole or, as explained in the 
next section, make use of the natural Brownian rotation 
of nanodiamonds in solution. 


In addition to the rotation of the externally applied 
magnetic field for nanodiamond powder, another option 
consists of the use of random Brownian rotations of nan¬ 
odiamonds in a solution. If the Brownian rotation is suf¬ 
ficiently slow to satisfy our adiabatic rotation condition, 
the polarized nuclear spins will maintain their polariza¬ 
tion in the direction of the externally applied magnetic 
field. The time scale of Brownian rotations is determined 
by 0 


Tb 


3 V H r] 
kT 5 


(34) 


where k is the Boltzmann constant, T is the temperature, 
Vh is the hydrodynamic volume of the particle and 7 
the viscosity of the surrounding carrier liquid. Note that 
the hydrodynamic volume is an effective volume that in¬ 
cludes both the true particle volume and the volume of 
a fluid that is displaced when the particle rotates due 
to particle-fluid interaction. A typical average hydro- 
dynamic diameter for nanodiamonds of 30 nm in water 
gives a Brownian relaxation time of tb — 9 /is, which 
is comparable with our polarizing cycle time. One can 
increase the Brownian relaxation time by using nanodi¬ 
amonds with a larger diameter, a more viscous fluid or 
lower temperatures. Additionally, surface coating can 
increase the average hydrodynamic diameter, i.e., after 
adsorption of ferritin the diameter of nanodiamond was 
determined to be 85 nm [4lj |. 

Here we assume that the hydrodynamic diameter of 
nanodiamond is 85 nm, and the Brownian relaxation time 
is tb — 205 /is. On average after a time tb, the nanodi- 
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amonds will have rotated to another random direction, 
and different NV centers will now participate in the po¬ 
larization dynamics of nearby nuclear spins while other 
NV centers will be inactive as they are effectively de¬ 
coupled from their neighboring nuclear spin due to the 
absence of a Hartmann-Hahn condition. When an NV 
center ” leaves” the coherent exchange range via a ran¬ 
dom rotation, the polarization of the polarized nuclear 
spins will still be transferred to the other nuclear spins 
by spin diffusion. Considering maximum of 20 degrees 
deviation, 6 % percents of the NV spins are involved at 
each given time, and each NV will be in this range after 
16 random rotations on average. 


IV. INTERNUCLEAR INTERACTION AND 
SPIN DIFFUSION 


So far, we have considered polarization of an individ¬ 
ual nuclear spin and have therefore ignored the effects 
of inter-nuclear coupling on both the polarization cycle 
and on the nuclear spin polarization diffusion towards 
more distant nuclear spins that are not directly inter¬ 
acting with the NV center. This section discusses these 
effects of nuclear dipolar interactions. Considering our 
off-resonant driving in the dressed state basis and the 
suitable rotating frame, the effective Hamiltonian of the 
whole system consisting of an NV spin and its nearby 
nuclear spins can be written as 


Htot ~ H nv + H h f + H d _ d + H Nu . (35) 


Here Hnv denotes the NV center effective energy Hamil¬ 
tonian, Hhf the effective electron-nuclear hyperfine inter¬ 
action, H d _ d the dipole-dipole coupling among nuclear 
spins and Hn u is the magnetic field splitting of nuclear 
spins 


Hnv 

Hhf 

H d - d 

Hnu 


2A(t)a x 

H - ^eff&z 5 

(36) 

5>(2« 

lx'. I?' 

X v 

+ 2 a z'Ji ) 

(37) 


V - 

- 3(1 i • e Vij )(Ij • e rij ) 

(38) 





YslnBlf, 


(39) 


with a x { and a z > denoting the hyperfine coupling strength 
between the ith nuclear spin and the electron spin, dij = 
j - 7n s n is the coupling strength between the i and j nu- 

77 r ij 

clear spins, with = |r^ | denoting the distance be¬ 
tween the two nuclear spin and e rij = . The i and 

j indices are summed over all nuclei. 

In order to estimate the effect of internuclear interac¬ 
tions on the polarization transfer from the NV center to 
the nuclei in the multiple nuclear spins case we consider 
5 nuclear spins and one electron spin, including hyper¬ 
fine and nuclear dipolar couplings. After each polariza¬ 
tion cycle, the electron spin is re-initialized in the state 


| — 1). The concatenated evolution of the nuclear spin 
density matrix is determined by an equation analogous 
to eq. ([33]) . The initial density matrix of the five unpo¬ 
larized nuclear spins can be written as po = I/2 5 , with I 
being a unit matrix of dimension 2 5 , and the time evo¬ 
lution is given by the operator U t = e ~ z J Htotdt ? w ith 
H to t denoting the total Hamiltonian eq. (l35|) . The po¬ 
larization is defined as in section II, P k = (I z r)/(I z r)o 
in which the subscript k presents different nuclear spins 
and P = \ Tl Pk • We assume in these simulations only 
nearest neighbour inter-nuclear spin coupling with rates 
d\2 = ^23 = 0(34 = c /45 = 2 kHz. 

The simulation results are presented in Fig. fTOh and 
b. We consider the case that the NV is strongly coupled 
to all five 13 C spins, a x / = (27r)0.7 MHz, a x ^ = (27r)0.5 
MHz, a x ' 3 = (27r)0.4 MHz, a x ^ = 0.32 MHz, and a x /. = 
(27r)0.2 MHz. The polarization is built up for all the 
five nuclear spins with different coupling to the NV spin, 
which matches our theory. In this case, the coupling be¬ 
tween the nuclear spins is too small to affect the evolution 
of the polarization dynamical built-up (Fig. H0h). Due 
to small interaction among the nuclear spins, there is al¬ 
most no difference between polarization in the presence 
or absence of the interactions among the nuclear spins, 
which means our scheme works in the frozen core. 

Now we would like to exemplify the benefits of inter¬ 
nuclear interactions which leads to nuclear spin diffusion 
and a considerable extension of the range of polariza¬ 
tion. Suppose a typical configuration of nuclear spins 
distributed around the NV. To estimate the polariza¬ 
tion efficiencies in a large system, we adopted the sim¬ 
ple spin temperature approximation, namely neglecting 
the nuclear spin coherence (42|, which provides a good 
estimation of polarization transfer efficiency when the 
sweep time step is sufficiently small. This approximation 
yields independent rate equations for each individual nu¬ 
clear spin and should give a conservative estimate for the 
achieved polarization. For the At sweep times described 
above, this approximation serves as a lower bound, as 
the achieved polarization in one cycle with our scheme 
is higher than the polarization achieved during the same 
time with small sweep steps. 

For the diffusion step, the interaction between the nu¬ 
clear spins is taken into account, where we use the Gaus¬ 
sian approximation for the nuclear spins. As is typical for 
paramagnetic centers, the NV center strongly affects the 
spin diffusion in nearby nuclear spins (the ” frozen core”) 
due to the energy mismatch caused by the S Z I Z hyperfine 
term. However, when the NV center is in the | m s = 0) 
state, this energy mismatch is suppressed and diffusion 
to external nuclear spins is allowed. To give a conserva¬ 
tive estimate, we assume that no diffusion in the frozen 
core takes place when the NV center is in the \m s = ±1) 
states. For this purpose, we simulate the NV center as a 
classical spin with the probability of being in each eigen¬ 
state given by the optical initialization corresponding to 
the NV orientation. 

As the nanodiamonds are in a solution, we assume 
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FIG. 10. The polarization buildup for five nuclear spins for Q e ff — 3.23 MHz. The nuclear spins are strongly interacting 
with the NV center spin, a x ^ = (27r)0.7 MHz, a x ^ = (27r)0.5 MHz, a x / = (27r)0.4 MHz, a x ^ = 0.32 MHz, and a x ^ = (27r)0.2 
MHz: a, The blue and red triangles denote the presence and absence of the interaction among the nuclear spins. Internuclear 
interactions do not lead to significant changes in the polarization dynamics because the dipole-dipole interaction of nuclear 
spins is much smaller than the electron-nuclear coupling, b, The polarization dynamics of the individual nuclear spins in the 
set of five spins represented by different markers within different colors, c, Numerical simulation of the polarization buildup 
in a large number of polarization cycles in a system consisting of an NV spin and around 548 13 C nuclear spins in a volume 
of 50000 lattice sites by using Brownian rotation. The internuclear dynamics is treated within the bosonic approximation and 
we simulate a polarization buildup time of 2 second. The system benefits from polarization diffusion between nuclear spins 
resulting in a high level of polarization. 


that the differences in polarization transfer for the dif¬ 
ferent orientations of the NV spins within the 20-degree 
deviation range can be neglected. Consider a random 
configuration of nuclear spins, with 548 nuclear spins in 
the vicinity of an NV center (lattice of 50,000 nuclei per 
NV center) whose orientation is within the 20-degree de¬ 
viation range. The net polarization is built up by us¬ 
ing Brownian rotations as shown in Fig. [T0b . There 
are two steps in the simulation. First, we assume that 
the orientation of the NV center satisfies 0 G [0°,20°] 
and we consider Q e ff = (27r)2.3 MHz, sweep rate of 
v = (27r)0.8 MHz //is and At = 70 /is. The polariza¬ 
tion transfer from the NV center to the neighboring nu¬ 
clear spins is achieved by using our combination of off- 
resonantly driven double quantum transition and the ISE 
technique. The polarizing cycles are repeated during the 
Brownian relaxation time tb = 205 fi s. Secondly, after 
the Brownian relaxation time, the nanodiamond is ran¬ 
domly rotated to another direction. As we discussed in 
the Sec. nm it takes roughly 3.5 ms for the nanodia¬ 
mond to rotate back to our maximal deviation range for 
polarization, and nuclear spin diffusion dominates dur¬ 
ing this time period (taking into account the NV effect 
on diffusion in the frozen core). In this step more distant 
nuclear spins can become polarized by diffusion. These 
two steps occur sequentially over a 2 second time range. 
In line with previous sections, the polarization is defined 
as = (I z t)/(I z t) o, where I z > = M is the number 

of the involved nuclear spins. Interestingly, the slope of 
polarization is close to linear. This demonstrates that in 
the long diffusion time, where the nanodiamond is not 
involved in the polarization cycles, we do not see a bot¬ 
tleneck caused by slow diffusion and significant polariza¬ 
tion is transferred from the NV spin to the nuclear spins 
within each step 1. The final net polarization reaches 


0.2, equivalent to about 110 13 C fully polarized spins are 
polarized within 2 second in our scheme. 


V. A LARGE ANGLE POLARIZATION 
SCHEME 


So far we have concentrated so far on a range of angles 
6 G [0°,20°] between NV center and the external mag¬ 
netic field. Note however, that this is not the only range 
in which efficient polarization transfer can be achieved. 
Indeed, we have another interesting range of the NV ori¬ 
entations, concentrated around the direction perpendic¬ 
ular to the magnetic field, i.e., 0 G [70°, 110°], which 
involves 34% the NV spins. Here the condition D(6) < 0 
leads to the Hamiltonian H _ in eq. m • Following steps 
that are analogous to those leading up to Hamiltonian eq. 
((201) we find for the subspace {|%_, 4-), |x+>t}} 


( U tt - a*' sin y \ 

ZJ _| W e// 2 2 | 

matnx - 'v -u ef f + *£ ) 


(40) 


Suppose now that the NV center is initialized in state 
| — 1) by using the off-resonant driving. Then we make 
use of the ISE technique to sweep adiabatically across 
the two possible Hartmann-Hahn resonance points of the 
dressed NV spin and nuclear spin pair. As a result the 
the population is transferred between |x_, \) and |x+, t), 
with the other states being unaffected. As before this 
allows the nuclear spins to be polarized but now in a di¬ 
rection opposite to the case 0 G [0°, 20°]. In the following 
we discuss the details of this case. 













14 


A. Initial polarization of the NV spins 

We begin with a discussion of the initialization of the 
NV center spins. As explained in previous sections op¬ 
tical pumping of the NV spins results in the \m s = 0) 
state in the NV frame, that is 

|0) e = cos0|O> + ^(e^| + 1) - 1>) 

where the states |0), | ± 1) are defined in the laboratory 
frame. For 0 E [70°, 110°] the population in state |0) is 
very small and the initial state is well approximated by 
A=(e l( ^\ + 1) — e~ l( ^\ — 1)) which, in the relevant subspace 
spanned by {| + 1), | — 1)}, is unpolarized. 

We address this problem by means of an adiabatic 
sweep. For concreteness, consider Q = (27r)20 MHz 
and an adiabatic sweep from A wm(U) = (27r)320 MHz 
to A wAf(tf) < — (27r)320 MHz. For a sweep range 
|co>i (tf) — LJi(ti)\ ~ (27r)640 MHz, the time for an adi¬ 
abatic sweep is of the order of 0.4 /is. Such a sweep in¬ 
duces an adiabatic population transfer from state | — 1) 
to state |0) so that the state of the NV center evolves to 

|0)@—> cos6*| - 1) + ^(e‘*| + 1) - e"^|0)) 
v 2 

where we ignore any dynamical and geometric phase that 
has been accumulated in the sweep. It is now crucial to 
note that the state |0) does not contribute to the polar¬ 
ization dynamics as it does not take part in the far de¬ 
tuned dynamics that is induced by the applied microwave 
fields. Hence the relevant quantity is the polarization in 
the subspace spanned by the states {| + 1),| — 1)}. Nor¬ 
malized by the population that is found in this subspace 
we find 



where N r = 1/ ( sm 2 ^ -hcos 2 9) is the normalized coefficient 
and P^ 2 ~ 0.85. 

B. The sweep range for ISE 

Owing to the different behaviour of the energy 
shifts D(9) and 5(9) we need to change the sweep 
range in the ISE. In order to achieve optimum po¬ 
larization transfer, we choose £2 = (27r)40.5 MHz, 

which results in a range of effective Rabi frequency 
Sleff e [(2 tt)2.3MHz, (2tt)3.6MHz\. The adia¬ 
batic slow passage needs to cross the range A E 
[—{2'k)2\MHz,(2'k)2\MHz]. For a sweep rate of v = 
(27r)6 MHz //is we estimate the required polarizing time 
as At > 7 /is. Apart from these parameter changes, the 
principles underlying our discussion of the polarization 
transfer in Section II are still valid. 


If we assume that the NV spin is initial polarized 
perfectly, we have P maa , > 0.2 in a coupling range 
a x > E [0.4,0.9] for \ft eff \ = (2tt)2.27 MHz in Fig. Efc, 
compared to a x > E [0.18,0.9] for \Q e ff\ = (27r)3.6 MHz 
in Fig. [6b. We can say that our off-resonant driving 
and ISE technique allow polarization transfer for all the 
NV orientations within 0 E [70°, 110°]. Consider the av¬ 
erage polarization of the NV spins in this case is about 
0.5 and twice effective Rabi frequency range here, apply¬ 
ing a singe sweep speed becomes not very efficient for 
all the NV orientations, especially for the nuclear spins 
with weak couplings to the NV spin. As we discussed, 
we have another choice, slow down the ISE sweep can 
involve the nuclear spins with weaker coupling to the NV 
spin, i.e., as shown in Fig. [8b, we can use another a speed 
rate v = (27r)0.8 MHz//is (At ~ 50 fi s) to make it pos¬ 
sible to polarized the 13 C spins with coupling strength 
a x > — (27r)0.1 MHz. 

In conclusion, off-resonant driving and the ISE tech¬ 
nique allow for polarization transfer between the NV spin 
and its surrounding nuclear spins, when the orientation 
of the NV spin matches 6 E [70°, 110°] but we need to 
sacrifice half the population to provide a high initial po¬ 
larization of the dressed NV spins. On the other hand 
the fact that around 34% of the NV centers participate 
in the polarization dynamics even without magnetic field 
rotation is an attractive feature. 


VI. DEPOLARIZING EFFECT 

Our approach to increase the fraction of NV center 
that participate in the polarization sequences by means 
of adiabatic rotations may potentially have an undesired 
side effect as nuclear spins that have been polarized for 
one magnetic field orientation may be depolarized again 
when transferred to a different magnetic field orienta¬ 
tion. Consider for example the case shown in Fig. [9] 
where, after polarization transfer is achieved for the blue 
spherical sectors, the magnetic field is rotated from in¬ 
direction to different direction Z '. If the NV spins in the 
blue sectors have not been polarized in advance and our 
combination of off-resonant driving and ISE achieves a 
Hartmann-Hahn condition at some instant, then the de¬ 
polarized state of the NV center will be transferred to 
nuclear spins, whose polarization reduces as a result. In 
the protocol for the case 0 E [0°, 20°] this depolarization 
is negligible due to the rapid increase of the detuning 
5(9) outside of this range. Only NV centers in the range 
0 E [20°,25°] may still experience undesired resonances. 
As these NV centers are still well initialized by means of 
optical pumping depolarization is negligible. 

The situation is different and indeed more complex for 
the second protocol that covers the range 0 E [70°, 90°]. 
NV centers with orientations in the range 0 E [35°,43°] 
are difficult to be polarized to the required state. Un¬ 
fortunately, they have almost the same energy distribu¬ 
tions \D(0)\ and 5(0) as the NV centers in the range 
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FIG. 11. The depolarization of the previously polarized nu¬ 
clear spin with B = 0.36 T, Q e ff — (27r)3.5 MHz, v = (27r)6 
MHz / /is and At = 10 ps. The coupling strength is given as 
a x t = (27r)0.3 MHz, (blue circles) and a x > — (27r)0.6 MHz 
(red triangles). 


6 G [70°, 110°], see Fig. [2j Our off-resonant driving and 
ISE technique scheme can lead to Hartmann-Hahn res¬ 
onances at some instant which then cause polarization 
exchanges between the unpolarized NV centers and the 
polarized nuclear spins. Thus, magnetic field rotations 
could lead to depolarization dynamics in the previously 
polarized nuclear spins for the protocol that is adapted 
to the case 0 G [70°, 110°]. 

We consider a simple system consisting of one NV and 
one nuclear spin as an example for studying the depo¬ 
larizing effect of this case. The depolarizing cycle time 
is as the same as our polarizing cycle time, since every 
initialization process is implemented on all the NV spins. 
Consider a nuclear spin that is initially polarized and in 
contact with an unpolarized electron spin p e o =1/2. The 
evolution of the nuclear spin density matrix is then given 

by 

P -> • • • Ultr e [Ul (I 1)4 I ® Peo)U?] ® PeoU? • • • , (42) 

/ s' 

matrix 1 j s time evolution operator. 

Hmatrix * s given by eq. (j20|) with the choice Q e ff = 
(27r)3.5 MHz and S' denoting that a magnetic field in in¬ 
direction is being considered. The polarization P of the 
nuclear spin is defined as in Sec. II. As shown in Fig. 
EU the nuclear spin is depolarized very rapidly. The de¬ 
polarization range O' G [35°,43°] includes 9% of the NV 
orientations while the range in which the polarization is 
active includes 34% of the NV centers. Hence, although 
adiabatic rotation of the magnetic field brings an unex¬ 
pected depolarizing effect, it is nevertheless quite attrac¬ 
tive as the fraction of NV spins that are contributing to 
the polarization transfer even without magnetic field ro¬ 
tation strongly outweighs the small fraction of 9% of the 
NV centers that suffer depolarization. 


VII. DISCUSSION 

To summarize, our polarization scheme for ensembles 
of nanodiamonds (see Fig. [7] and Fig. [9]): 

Step 1: A moderate to strong magnetic field is applied 
along Z-direction (of order B = 0.36T). A polarization 
cycle consists of optically pumping and initialization of 
the NV center spins, and the subsequent application 
of off-resonant microwave driving and the integrated 
solid effect to realise a robust adiabatic sweep across 
Hartmann-Hahn resonance points. This achieves the 
coherent polarization transfer from the NV center spins 
to the nuclear spins with a range of orientations of 
the NV center relative to the magnetic field. This 
polarization cycle is repeated several times. 

Step 2: Adiabatic rotations (of the magnetic field or of 
the nanodiamonds) change the magnetic field orientation 
relative to the NV centers to a new Z'-direction, while 
the nuclear spins follow the direction of the magnetic 
field. 

Step 3: Repeat step 1 and step 2. 

Apart from control imperfections, the achievable 13 C- 
bath spin polarization will be limited by the spin-lattice 
relaxation time T\ of the nuclear spins. In recent experi¬ 
ments, by using electron paramagnetic resonance (EPR), 
solid-state NMR, and DNP techniques, a variety of nan¬ 
odiamond samples, of varying manufacturing processes 
and particle sizes, were examined (l8j |. Several types of 
nanodiamonds were identified that exhibit very long nu¬ 
clear spin lifetimes of several minutes. Additionally, this 
polarization transfer rests on NV spins having a suffi¬ 
ciently long rotating-frame spin-lattice relaxation time 
Tip in nanodiamonds during the adiabatic passages that 
are required for the implementation of ISE. Using a 
strong MW field driving, T\ p can be estimate as 100 
/is [39|, |43|, which is much longer than our sweep time 
in one polarization cycle. Assuming a polarization cycle 
duration of At =10 /is, when rotating the magnetic field, 
17 rotations of the magnetic field are required to involve 
all the NV spins. For a nuclear spin lifetime of more than 
2 minutes, this would allow the execution of around 10 5 
polarization cycles for each NV spin involved before the 
nuclear spins start to relax. For nanodiamonds in a so¬ 
lution, each nanodiamond is involved in the polarization 
cycles roughly once every 17 random Brownian rotations. 
Assuming tb — 205 /is, with a pure nuclear diffusion pe¬ 
riod of roughly 3ms after every Brownian relaxation time 
tb , 2 minutes lifetime of the nuclear spin would allow in 
excess of 10 5 polarization cycles for every NV center. 

We can now obtain a rough estimate of the total 
amount of the net polarization that could be achieved un¬ 
der realistic experimental conditions. Consider a 1 mm 3 
nanodiamonds powder, dissolved in a solution. The nan¬ 
odiamonds are assumed to have a diameter of 85 nm, 
with an NV center concentration of C = 2 x 10 18 cm -3 
(see [44] and references therein for achievable NV con¬ 
centrations). In the limit of a uniform spatial distribu- 
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tion, this gives about 642 NV spins per nanodiamond, 
with the total amount of NV spins given by 8.35 x 10 15 
NV spins. Consider natural abundance of the 13 C spins 
( 1 . 1 %), the total nuclear spins in this powder of nanodi¬ 
amond is approximately 8.18 x 10 18 [43]. Note that 13 C 
enriched nanodiamonds will enable polarization of more 
nuclear spins with each NV center and hence a larger 
overall level of polarization. Taking account of the high 
degree of polarization, in excess of 0.3 as demonstrated 
in Sec. □30 we find that a total amount of polarization 
equivalent to 1.6 x 10 18 nuclear spins should be achievable 
in such a sample within 2 seconds. 


VIII. CONCLUSION 


In conclusion, in the high magnetic field limit, we pro¬ 
pose a scheme achieving macroscopic levels of 13 C nu¬ 
clear spin polarization in randomly oriented ensembles 
of nanodiamonds, realised as powder or solutions. To 
address the lack of a common quantization axis for the 
NV centers, an off-resonant microwave drive realises a 
resonant double quantum transition and this together 
with the integrated solid effect enable microwave cou¬ 
pling and control of NV spins whose orientation deviates 
from the external magnetic field by less than 20 degrees. 
Matching the effective Rabi frequency of the off-resonant 
double quantum transition of the NV centers to the nu¬ 
clear spin Larmor frequency, enables near resonant cou¬ 
pling between NV electron spins and the nuclear spins, 
and the transfer of polarization from the initialized NV 
spins to the neighbouring 13 C nuclear spins. Addition¬ 
ally, the effect of nuclear dipole-dipole interactions allows 
the weakly coupled nuclear spins to be polarized by spin 
diffusion. Adiabatic rotation (of the magnetic field or 
sample) can then extend the polarizing scheme to more 
NV spins, achieving a net polarization of all nuclear spins 
via our scheme. By using our polarization scheme, room- 
temperature, optical based polarization of the nuclear 
spins is made possible. 

These results introduce several exciting opportunities 
for the application hyperpolarized nanodiamonds, espe¬ 
cially in the biomedical sciences. Diamond nanoparticles 
are biocompatible, exhibit no in vivo toxicity (46j, and 
can be attached to a wide range of specific proteins, pep¬ 
tides or antibodies 0 while maintaining long T1 times. 
Thus, hyperpolarized nanodiamonds present an exciting 
platform as MRI probes for molecular imaging. 

I 

/ D(6) + 7 e B z 

-Gt 

V 


APPENDIX: EFFECTIVE HAMILTONIAN OF 
AN NV CENTER AT THE HIGH MAGNETIC 
FIELD LIMIT 

In our work we assume that in the laboratory frame 
of reference a strong magnetic field is applied along the 
^-direction, 7 e B z D. In the following we derive the 
representation of the NV center Hamiltonian eq. [3] in the 

^0 



FIG. Al. Two related rotations of the coordinate system are 
used to map from the natural orientation zq - axis of the NV 
spin to z-axis which is defined by the high magnetic field. 

laboratory frame of coordinates. To this end we need to 
transform the zero-field splitting tensor, which is diago¬ 
nal in the natural coordinate system defined by the NV 
center axes, to the laboratory frame of reference. We as¬ 
sume that the angle between the magnetic field, i.e. the 
^-direction, and the NV-axis, the ^-direction is given 
by 6 and that the angle between the laboratory frame y- 
direction and the NV center ye directions is given by in 
the x-y plane, as shown in Fig. IA11 The spin operators 
in the NV system of coordinates, S Xe , S ye and S Ze , are 
then related to the spin operators S x , S y and S z in the 
laboratory coordinates by 

S X0 = sin 0S Z + cosQ(cos(l)S x — sin (f>S y ), 

S ye = cos (j)S y + sin (j)S x , 

S Ze = cos 0S Z — sin # (cos — sin 4>S y ). 

The Hamiltonian of the electron spin 

H e N = (—£> + E)s xe • S x „ + (-Id - E)S yt • Sy, 
+ ( 2 -D)S Ze -S Ze + le BS z , 

can then be rewritten in the basis defined by the eigen¬ 
states of the S z operator of the electron spin, {I + 
1), |0), I - 1)} as 


-Gi G 2 

0 G\ 

G* D{0)- le B z 


TTG _ 

n NV ~ 
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with 

_ D{ 1 + 3 cos(2#)) + 3E(1 - cos(2 0)) 

= -4-’ 

(ZJ — E 1 ) sin # cos Oe 1 ^ 

■ = — ?! — ■ 

^ [D + 3E + (E - D) cos 26>]e 2 3 4 5 ^ 

G 2 =-4-• 

There exist two uncontrollable components, the angles 
6 and <f> as well as the local strain E, which modifies the 

I 


spin level structure. 

As we know that a high magnetic field suppresses the 
effect of the off-diagonal terms in H'^ v . Assume that the 
Hamiltonian is of the form H r ^ v = H'nyo + eV, where 
is the diagonal part of the Hamiltonian and eV 
the off-diagonal part. In our case the magnetic field is 
sufficiently large for us to assume that the off-diagonal 
part is a weak perturbation. By using the Schrieffer- 
Wolff transformation in condensed matter, the second 
order corrections due to the off-diagonal terms eV can 
be obtained as 


WKi\f3) 


7 W#f) 

2 E a — Ei 


Ei-Ep ^ 


in which H^ is defined as the corrected Hamiltonian. By 
simple calculation, we find 


(+!|#mI + 1) — 

(-M h m\ ~ 1 ) 


_ , |g 2 | 2 

7 e B + D(6) 2 7e H’ 


\Gi 


\Gi 


\G2 


- le B + D(0) 2 le B' 

After shifting the zero of energy, the effective Hamilto¬ 
nian in the laboratory frame is given by 


K ff = 7 eBS z 


D(0)S 2 Z +5(0)S Z , 


in which 


m = 


7 e^|Gi | 2 |G 2 | 2 

( 7e H) 2 - [D(0)]> ^ 2 le B' 


Therefore a strong magnetic field suppresses the first or¬ 
der effect of the off-diagonal terms in H^ v but results 
in a second order modification of the diagonal elements 
5(6). This is the Hamitonian as in eq. (j3j). 


In order to estimate the validity of this approximation 
in the high magnetic field limit we have prepared the sys¬ 
tem in state \m s = 0) and observed the subsequent time 
evolution of its population under Hamiltonian H ^ in Fig. 
IA2l (see caption for simulation parameters). It is evident, 
that the population of the initial state m s = 0 remains 
essentially constant which implies in turn that the effect 
of the off-diagonal elements in Hamiltonian H r ^ v are in¬ 
deed negligible for the moderate to high magnetic field 
case that is relevant to our work. 
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